Entanglement in a Valence-Bond-Solid State 
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We study entanglement in Valence-Bond-Solid state, which describes the ground state of an AKLT 
quantum spin chain, consisting of bulk spin-l's and two spin-l/2's at the ends. We characterize 
entanglement between various subsystems of the ground state by mostly calculating the entropy 
of one of the subsystems; when appropriate, we evaluate concurrences as well. We show that the 
reduced density matrix of a continuous block of bulk spins is independent of the size of the chain 
and the location of the block relative to the ends. Moreover, we show that the entanglement of the 
block with the rest of the sites approaches a constant value exponentially fast, as the size of the 
block increases. We also calculate the entanglement of (i) any two bulk spins with the rest, and 
(ii) the end spin-l/2's (together and separately) with the rest of the ground state. For example, 
we show that (i) any two bulk spins become maximally entangled with the rest of the ground state 
exponentially fast in their separation distance, (ii) the two end spin-l/2's share no entanglement, 
and (iii) each end spin-1/2 is maximally entangled with the rest. 

PACS numbers: 75.10.Pq, 03.67.Mn, 03.65.Ud, 03.67.-a 
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There is considerable current interest in quantifying 
entanglement in various quantum systems. Entangle- 
ment in spin chains, correlated electrons, interacting 
bosons and other models was studied B, B , B , B, B, B , 
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|2J, |25|, |26L |27| . Entanglement is a fundamental measure 
of how much quantum effects we can observe and use, 
and it is the primary resource in quantum computation 
and quantum information processing [ij, ^ . Also entan- 
fement plays a role in the quantum phase transitions 
I , and it has been experimentally demonstrated that 
the entanglement may affect macroscopic properties of 
solids jailal- 

In this Letter, we will study a spin chain introduced 
by AfHeck, Kennedy, Lieb, and Tasaki (AKLT model) 
I2M I23 . The ground state of the model is a unique pure 
state. It is known as Valence-Bond Solid (VBS), and 
plays a central role in condensed matter physics. Hal- 
dane [3l| conjectured that an anti-ferromagnetic Hamil- 
tonian describing half-odd-integer spins is gap-less, but 
for integer spins it has a gap. AKLT model describing 
interaction of spin-l's in the bulk agrees with the con- 
jecture. An implementation of AKLT in optical lattices 
was proposed recently |32| , and the use of AKLT model 
for universal quantum computation was discussed in [SOJ . 
VBS is also closely related to Laughlin ansatz [S^] and 
to fractional quantum Hall effect [3J|. 

We investigate the seminal AKLT model from the new 
perspective of quantum information, and evaluate the en- 
tanglement (in terms of entropy) of various subsystems 
of the VBS. The results and the methodologies adopted 
herein have several implications from the perspective of 
both quantum information and condensed matter. For 
example, while the entanglement in spin chains with 
periodic boundary conditions has been studied exten- 
sively, our results provide entanglement calculations for 



spin chains with open boundary conditions. For criti- 
cal gap-less models conformal field theory describes the 
entanglement l2^|2^|23. For gapped models G. Vidal, 
J.I.Latorre, E.Rico and A.Kitaev Q conjectured that the 
entropy of a large block of spins reach saturation. We 
confirm this for the AKLT model, and find that the en- 
tropy of a large block of bulk spins is close to two. This 
means that the block can be in four different states, and 
hence, the Hilbert space of states of the large block of 
bulk spins is four-dimensional. Our results also show 
that the entanglement correlation of VBS state is short- 
ranged, which provides a good understanding why the 
density matrix renormalization group (DMRG) method 
[4J| works so efficiently for VBS states; see 45] for recent 
developments. 

The AKLT model consists of a linear chain of N spin- 
l's in the bulk, and two spin-l/2's on the boundary. We 
shall denote by Sk the vector of spin-1 operators and by 
Sft spin-1/2 operators, where 6 = 0, TV -I- 1. The Hamilto- 
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1 is: 




N-l 


H - 


k=l 



SkSk+1 + -{SkSk+i) + ttq,! + 7rjv,jv+i- (1) 



The boundary terms tt describe interaction of a spin 1/2 
and spin 1. Each term is a projector on a state with spin 

3/2: 

Ti'o,! =0(1 + soSij , 7rjv_Ar+i =0(1 + ^n+iSn) ■ (2) 

The ground state of this model is unique and can be 
represented as J2^ |23 : 

|G)(®LiPm.)I*'>oiI*~>12 • • • l*->^^+i- (3) 

Here P projects a state of two qubits on a symmetric 
subspace, which describes spin 1. In the formula above 



I^~)(l Ti) — I iT))/^^ represents a singlet state, and the 
subscripts represent the two parties the singlet is shared 
between. We have tried to keep our notations as close to 
those in the paper jl3|- We can use the following figure 
to visualize the ground state: 



I*-) I*-) I*-) ll*-) ll*-) I*-) I*- 
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A black dot represents spin--l, and spin-l's are denoted 
by circles. To begin with, each bulk site, k (where 
1 < fc < N) shares one singlet state |^~) (represented 
by a line) with its left and right neighbors. Thus at 
each bulk site, k, we start with two spin-l/2's labeled by 
{k,k) and then the spin-l's are prepared by projecting 
the two spin-l/2's (4-dimensional space) on a symmet- 
ric three dimensional subspace of spin 1 (3-dimcnsional) . 
The system has open boundary conditions, and the two 
ends are numbered as sites (before projection, this site 
shared a singlet with site 1) and A^ -I- 1. 

There is an upper bound on the entropy of a block 
of L spins. Before projection, the entropy is equal to 
2, since the boundary intersects two singlet states. Since 
the local projections will only decrease the entanglement, 
we expect that the entropy of a block of L spins to have 
an upper bound of 2. 

In order to calculate the reduced density matrices of 
various subsystems of the ground state \G) (see Eq.Ol), it 
is more convenient to express it in a different form based 
on the singlet chain shown in the preceding figure and 
the figure below. 
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Let us first consider a chain of two singlet states, \'9~)ab 
and \'^~)bc'- -^ i^ '^^ ^ite #1, (B, B) is in site #2, and C 
is in site #3. The combined state can then be expressed 
as follows: 



1 ^ 



a=0 



(8)/a (8) (CTa)c) I* )sbI* )aC, 



(4) 



where both / and ctq represent the identity operator, 
o'It<J2, cr^i are the Pauli matrices, and '*' means complex 
conjugation. By entanglement swapping similar to tele- 
portation [33, party #2 can perform a Bell state mea- 
surement on [B, B), and then communicate the results 
of measurements to party ^1 or #3. Then one of them 
can perform a unitary transformation locally, and finally 
a maximally entangled state will be shared by them. A 



multi-dimensional generalization of this can be found, for 
example in [36^. 

Eq.l^J can be generalized to a chain of singlet states. 
First, define quantum states |a) = (-l)^+"(/(g)(T*)|*~). 
Thus, |0) is the singlet state with spin 0, while other three 
states |1), |2), |3) form the symmetric subspace of spin-1 
(within a phase). Repeatedly using the relation I^J and 
with the help of the property presented later in the proof 
of our theorem, we obtain: 
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\ai) 



ai ,---,aiv— 



|aA')(4® (CTojv ••■Crai)w+l)|* )o 



0,iV-(-l- 



(5) 



The quantum states \ai) are orthonormal states at lattice 
site {i,i). Thus, by projecting the quantum state on the 
symmetric subspace spanned by the states |1), |2), and 
1 3), the ground state of AKLT model can be rewritten as 
12, 371: 



\G) 



1 



3JV/2 



E 



\ai)--- 



ai ,---,ajv— 1 
|aAf)(4«'(0'a„ •••CTai)N+l)|*")o,Ar+l. (6) 



It follows directly from Eq.® that the reduced den- 
sity matrix of spin-1 at any bulk site k (recall that 
k^ l,...,N) is: 



Pi 



1 ^ 

Tri,...{fe}...,Ar^o,JV+i|G')(G| = - E l"fe>("fc|, (7) 



oik=l 



where the trace is taken over all sites (including the two 
ends), except site number k. We see that all one-site 
reduced density operators in the bulk are the same: the 
identity or the maximally-disordered state in the spin-1 
space. Thus, the single-site reduced density matrices are 
independent of the total size of the spin chain A^, and of 
the distance from the ends (i.e., k or N — k). For the 
more general case, we have the following result: 
Theorem: Consider the reduced density matrix of a con- 
tinuous block of spins of length L (not including the two 
boundary 1/2-spins), starting from site k and stretching 
up to k -\- L — 1, where fc > 1 and k -\- L ~ \ < N (thus, 
1 < L < N) in the VBS ground state ^. Then, all these 
density operators are the same, and independent of both k 
(i.e., the location of the block) and of N (the total length 
of the chain). Thus, the reduced density matrix depends 
only on L, the length of the block under consideration. 

The proof is based on the following relations: Define 
|$+) = (I TT) + I ii))/V2, we know that |$+) = {-i){a2® 
/)|^^). For a unitary operator U, we have the property 
(C/® [/*)!$+) = |$+). Then (C/i C/2)|$+) = (f/iC/| "^ 
/))|(()+) = (/(g)C/2C/f)|<l>+), where Ui,U2 are two unitary 
operators (the super-index t denotes the transposition). 



By using these relations, we can prove that: 

Tro,Ar+i(/®t/ll^C/2)|*">(*"|(/®C/l^'C/2)^ 

= Ti^N+iil ®V)\<^+){^+\{I (g>Vy . (8) 

By repeated apphcations of this relation, and considering 
the ground state © , the reduced density operator of any 
continuous block of spins of length L is 



PL 



l-^K)(al|...K)KI 



3^ 

a,Q' 

TTo,N+i{I(^V)\^+){^+\{I(g>Vy. (9) 

where V = aaj^ ■ ■ ■ Ua^ , V = aa' • ■ ■ Cq,' . This operator 
only depends on L. This completes our proof. 

Our aim is to calculate the entanglement of the VBS 
state. For a pure bi-partite state \'iP)ab, the entangle- 
ment between spatially separated parties A and B is 
S{pa) = S{pb), where PA(s)TrB(A)|i/')(f/'| are the re- 
duced density operators and S{p) — — Trplogp is the von 
Neumann entropy, where we take the logarithms in the 
base 2. For example, it follows from Eq. {Tj) that the en- 
tropy of the one-site reduced density operator in the bulk 
is S {pi{k)) — log 3. This entropy describes the entangle- 
ment between site number k in the bulk (considered as 
one party) and the rest of the ground state (considered 
as the other party). The space of spin-1 is three dimen- 
sional, so log 3 is the maximum of the entropy. So we 
proved that in the VBS state (jHl, each individual spin 
in the bulk is maximally entangled with the rest of the 
ground state. Later in the paper, we shall see that this 
is also true for the boundary spin-l/2's. 

Since the reduced density operator of a continuous 
block of L spins is independent of the total size, N , of 
the spin chain, we can consider the case where L = N , 
i.e., we consider a chain of L spin-l's with one spin-1/2 
at each end. Now the reduced density operator of two 
end spin-l/2's takes the following form: 



1 



3 



P'L = ^ 2^ (/(»a„^---a„J|* )(* |x 

Ql,---,QL = 1 



t _ 



x(/(g)crQ^ ■■■aai) 
= \{l~p{L))-I + p{L)\^-){^-\. 



(10) 



Here p{L) — (—1/3)^ and / is the identity in 4 dimen- 
sions. Since the ground state © is pure, the entropy of 
the block of L bulk spin-l's is equal to the entropy of the 
two ends. So we have 

Sl = S{pl) = S{pi) = 
o , 3(l-p(i)) 



1 + ML) 



log(l-p(L)) 
log(l + 3p(L)). 



(11) 



As expected, Sl <'2 and approaches two 2 exponentially 
fast in L: S'l ~ 2 — (3/2)p(L). This is also clear from 
(110(1 : the reduced density operator approaches the iden- 
tity in the 4-dimensions exponentially fast. Consider the 
numbers: 



Sl = 1.58496 S2 = 1.97494 ^3 
5*4 = 1.99969 55 = 1.99996 Sq 



1.99695 



(12) 



Note that the correlation function of local spins decays 
equally fast: 



<5l5i>~(-1/3)^=p(L), 



(13) 



see [291,134 

Next we shall study the entropy of tv^ro spin-l's sep- 
arated by M sites in the bulk. That is we calculate the 
entanglement between two two bulk spin-l's and the rest 
of the spin-l's and the two spin-l/2's. We still can show 
that the reduced density operator does not depend on 
the total size of the chain, N, and prove that: 



P2iM) = -il-p{M))I+piM)p2, 



(14) 



where p{M) = (—1/3)^^ and p2 is the two-site reduced 
density operator of nearest neighbors, i.e. the case M — 
0, and the operator / is the identity in nine-dimensions. 
The nearest neighbor two-site reduced density operator 
can be written explicitly: 



P2 



= J[E l«)(/3|®l«>(/3| 



a,/3=l 



+ J2i\a){a\ |/3)(/3| - |a)(/3| ® |/3)(a|)]. (15) 

So we can calculate the entropy of two spins at distance 
M: 



S^iM) = 21og3 - -(1 -p(Af))log(l -p(M)) 



{l+p{M))log{l+p{M)) 



--(H-2p(A/))log(l + 2p(M)). 



(16) 



We see that S2{M) also approaches the maximum value 
(since the dimension is 9, the maximum entropy is 2 log 3) 
with the exponential rate defined by local correlations 
(HSJ. Note that S2 = 5'2(0) (seeEq.lHIl) and (O). How- 
ever, for M > 1, S2{M) quickly exceeds 5*^. We also can 
calculate the concurrence (another measure of entangle- 
ment [33). We shall use the generalized concurrence in 
higher dimensions J39j . Two concurrences corresponding 

to 5"^ and S2{M) are equal to Cl = I - P^{L) = I - -j- 



and C2{M) 



1 



1 



p\M) 



1 



1 



9^ 
r, respectively. 



6" '"'' 6-9*^' 

They look similar because the entanglement of the block 



also represents the entanglement of two ends with L bulk 
spins. 

Now we turn to the analysis of entanglement of 
boundary spins. We start from the reduced density op- 
erator of one boundary spin. We can prove that it is the 
identity matrix in two-dimensions. This shows that the 
end spin-l/2's are maximally entangled with the rest of 
the ground state, and has an entropy of 1. 

The density operator of two ends pj^ (^^^ Ea. ltTUll ') de- 
pends on the total size of the lattice N . In the simplest 
case, pi = [I — |\l/~)(^~|)/3. This is a separable state. 
Actually it is separable for any N . So, there is no en- 
tanglement between the two ends. If the size of the spin 
chain N increases, p^ approaches quickly the identity 
matrix in four dimensions. In Eo.f lll|) . replacing L by 
N in Sl we get the entanglement between the two ends 
(one subsystem) and all N bulk spins (another subsys- 
tem). It means that two ends considered as a subsystem, 
are maximally entangled with the bulk spins if N is large. 

Next we consider two-site reduced density operator 
with one spin in the bulk and another spin at an 
one end. It is enough to put the end spin at site 0, 
and the bulk spin at the site [M -1-1) (the range is 
M = 0,...,N — 1). We can calculate the reduced den- 
sity operator as above: 



I' 



p{M) 



P2(0,Af + l, g. , g 

+ |2)(l|®z(|l)(l|-|0)(0|) 

+ |3)(1|®(|0)(1|-|1)(0|) 

+ |3)(2|®(-z)(|0)(l| + |l)(0|)]. 



|l>(2|®z(|0)(0|-|l)(l|) 

|1)(3|«(|1)(0|-|0)(1|) 
|2)(3|®z(|0)(l| + |l)(0|) 
(17) 



where p{M) = (—1/3)*^. First we consider if this state is 
separable. Since it is the 2 x 3-dimensional case, we can 
use Peres- Horodecki criterion ,4(i,,41j. We find that when 
M ~ 0, the state is entangled. For M 7^ it is a separa- 
ble state. So, we know that the end spin 5 is entangled 
only with its nearest neighbor (spin-1). Secondly, we can 
study the entropy of this state, it is: 



1 



log 6 



(l-p(A/))log(l-p(Af)) 



(l + 2p(M))log(l + 2p(Af)). 



Similar to other entropies presented above, it approaches 
the upper bound log 6 with the same exponential speed, 
defined by local correlations (|13|l . The concurrence cor- 
responding to this entanglement is: C{0,M + 1) = 



In summary, the entanglement properties of the VBS 
state can be listed as follows: (1) each individual spin is 
maximally entangled with the rest; (2) the entanglement 
of a block of spins of length L with the rest gets to a 
constant value exponentially fast with L; (3) the entan- 
glement of any two bulk spins gets maximal exponentially 
fast in their distance; (4) each individual boundary spin 



is maximally entangled with its nearest neighbor and not 
with the other bulk spins and the other boundary spin; 
and (5) each individual boundary spin and another in- 
dividual bulk spin are entangled with the rest, and the 
entanglement gets maximal exponentially fast with the 
distance between the boundary spin and its bulk part- 
ner. 

In the future it will be interesting to calculate the 
entropy of a subsystem of 2 and 3-dimensional AKLT 
model; in fact, we are planning to study the entangle- 
ment for AKLT on arbitrary graphs using the results of 
[431 . We believe that it will be useful for universal quan- 
tum computation, as in |43| . 
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